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ABSTRACT 


This report presents the development of the Gain Weighted Eigenspace Assignment 
methodology. This provides a designer with a systematic methodology for trading off 
eigenvector placement versus gain magnitudes, while still maintaining desired closed-loop 
eigenvalue locations. This is accomplished by forming a cost function composed of a scalar 
measure of error between desired and achievable eigenvectors and a scalar measure of gain 
magnitude, determining analytical expressions for the gradients, and solving for the optimal 
solution by numerical iteration. For this development the scalar measure of gain magnitude 
is chosen to be a weighted sum of the squares of all the individual elements of the feedback 
gain matrix. An example is presented to demonstrate the method. In this example, 
solutions yielding achievable eigenvectors close to the desired eigenvectors are obtained 
with significant reductions in gain magnitude compared to a solution obtained using a 
previously developed eigenspace (eigenstructure) assignment method. 


1.0 INTRODUCTION 

The Direct Eigenspace Assignment (DEA) method (Davidson 1986) is currently being 
used to design lateral-directional control laws for NASA's High Angle-of-Attack Research 
Vehicle (Davidson 1992). This method allows designers to shape the closed-loop response 
by choice of desired eigenvalues and eigenvectors. During this design effort DEA has been 
demonstrated to be a useful technique for aircraft control design. The control laws 
developed using DEA have demonstrated good performance, robustness, and flying 
qualities during piloted simulation. These control laws are scheduled for flight test at 
NASA Dryden Flight Research Center in 1994. 

During the control law design effort, one limitation of the DEA method became 
apparent. Using DEA the designer has no direct control over augmentation gain 
magnitudes. Often it is not clear how to adjust the desired eigenspace in order to reduce 
individual undesirable gain magnitudes. To reduce undesirable gain magnitudes the 
designer must rely upon a strong physical insight into the dynamics or is forced to iterate 
on the design. 

This report presents the development of an eigenspace (eigenstructure) assignment 
method that overcomes this limitation. This method, referred to as Gain Weighted 
Eigenspace Assignment (GWEA), allows a designer to place eigenvalues at desired 
locations and trade-off the achievement of desired eigenvectors versus feedback gain 
magnitudes. 

This report is organized into four sections. Background information on how 
eigenvalues and eigenvectors influence a system's dynamic response and a review of the 
Direct Eigenspace Assignment methodology is presented in the following section. The 
development of the Gain Weighted Eigenspace Assignment methodology is presented in 
section 3. Concluding remarks are given in the final section. 



2.0 BACKGROUND 


This section presents a review of how eigenvalues and eigenvectors influence a system's 
dynamic response and a review of the Direct Eigenspace Assignment methodology. 


Eigenvalues, Eigenvectors, and System Dynamic Response 

The eigenvalues and eigenvectors of a system arc related to its dynamic response in the 
following way. Given the observable and controllable linear time-invariant system 

x = Ax + Bu 

and output equation 

y = Cx 

where x e R", u e R m , and y e R* . 

The Laplace transform of equation (2.1a) is given by 

sx(s) - x(0) = Ajc(s) + Bu(s ) 

x(s) = [si - A]~' x(0) + [si - A]" 1 Bu(s) (2.2b) 

Solution of equation (2.1a) is given by taking the inverse Laplace Transform of equation 

x(t) = .£ -l {[s/ - A]''}x(0) + £~'{[sl - A] -1 Rw(j)} (2.3) 

Noting that 

£~'{[sl - AY 1 } = e Al (2.4) 

the solution of (2.3) is (Brogan 1974) 

t 

x(t) = e A, x(0) + J e M, - r) Bu(T)dr (2.5) 

0 

and system outputs are 


(2.1a) 


(2.1b) 


(2.2a) 


/ 

y(t) = Ce A ‘ JC(0) + J Ce M, - x) Bu(r)dr 


( 2 . 6 ) 



The system dynamic matrix, A , can be represented by 


A = VAV l = VAL (2.7) 

where Vis a matrix of system eigenvectors, L is the inverse eigenvector matrix, and A is a 
diagonal matrix of system eigenvalues. Given this result, e^'can be expressed by 

e Al = Ve M L = X v, e X,t l. j (2.8) 

>=i 

where Xj is the j^ 1 system eigenvalue, is the j^ 1 column of V ( j^ 1 eigenvector of A ), and lj 
is the jth row of L ( j^ 1 left eigenvector of A ). Equation (2.6) can then be expressed as 

y(t) =XCv/' (,) / y x(0) + X Cvj fe x ' < " ti l J Bu( r) dr (2.9) 

y=l j=\ 0 


Noting that 


Bu(t) = ^b k u k (t) (2.10) 

Jt = l 

where b k is the k™ column of B and u k is the k 1 * 1 system input, the system outputs due to 
initial conditions and input u k is given by 


y(t) =£Cv y e l ’ u) lj x(0) + X X Cv j 1 A \e l>(, ~ r) u k (r) dr (2.1 1) 

j = 1 j=\ *=1 0 


The i 1 * 1 system output is given by 


?,(<>=£ 

y= i 


c i v j e 


A,(») 


^(0) + X X WA ]e l ' (, - T) u k (r)dr (2.12) 


;=1 k =\ 


where c ( is the i^ 1 row of C. In the case of initial conditions equal to zero, the i^ 1 output is 
given by ' 


*<') = XE R» k )e X,( '- x) u k {r)dr (2.13) 

y=i *=i o 

where R~ k = c- lj b k . In this expression, R^ k is the modal residue for output i, associated 
with eigenvalue j, and due to input k. 

Given an impulsive input in the k^ 1 input, equation (2.13) reduces to 
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XW = XX^« J ' < ' ) (2.14) 

i=\ *=i 

As this expression shows, for an impulsive input, a system's dynamics are dependent on 
both its eigenvalues and its eigenvectors. The eigenvalues determine the natural frequency 
pd damping of each mode. The eigenvectors determine the residues. The residues are an 
indicator of how much each mode of the system contributes to a given output. 


Direct Eigenspace Assignment Methodology 

One possible approach to the aircraft control synthesis problem would be to synthesize 
a control system that would control both the eigenvalue locations and the residue 
magnitudes associated with undesirable modes in certain outputs. Since the residues are a 
function of the system’s eigenvectors this naturally leads to a control synthesis technique 
that involves achieving some desired eigenspace in the closed-loop system ( Moore 1976; 
Srinathkumar 1978; Cunningham 1980; Andry 1983; Smith 1990). Davidson and Schmidt 
(Davidson 1986) explored this approach by using Direct Eigenspace Assignment (DEA) to 
synthesize flight control systems for flexible aircraft. DEA is a control synthesis technique 
for directly determining measurement feedback control gains that will yield an achievable 
eigenspace in the closed-loop system. For a system that is observable and controllable and 
has n states, m controls, and / measurements; DEA will determine a gain matrix that will 
place / eigenvalues to desired locations and m elements of their associated eigenvectors to 
desired values + . If it is desired to place more than m elements of the associated l 
eigenvectors, DEA yields eigenvectors in the closed-loop system that are as close as 
possible in a least squares sense to desired eigenvectors. The following section will present 
the development of the DEA synthesis technique. A more detailed development can be 
found in Davidson 1986. 


Direct Eigenspace A ssignment Formulation 
Given the observable, controllable system 

x = Ax+Bu (2.15a) 

where x e R" and u e R m , with system measurements given by 

z = Mx + Nu (2.15b) 

where z e R^ 

The total control input is the sum of the augmentation input u c and pilot's input u p 

« = « P +« C (2.16) 


t This assumes / > m. For a general statement and proof of this property the reader is referred to 
Srinathkumar 1978. 



The measurement feedback control law is 


u c =Gz (2.17) 

Solving for u as a function of the system states and pilot's input yields 

« = [/„- GN]-' GMx + [I m - GN]' u p (2. 18) 

The system augmented with the control law is given by 

x = (A + B[I m - GNY' GM )x + B[I m - GN]~' u p (2. 19) 

The spectral decomposition of the closed-loop system is given by 

(A + B[I m - GNY 1 GM )v, = A,v, (2.20) 

for i = 1 n where A, is the i” 1 system eigenvalue and v, is the associated i^ 1 system 

eigenvector. Let w, be defined by 


w, = [/„ - GN]~' GM v, (2.21) 

Substituting this result into equation (2.20) and solving for v,- one obtains 


v ( .=[A 1 / B -Ar 1 flw, (2.22) 

This equation describes the achievable i** 1 eigenvector of the closed-loop system as a 
function of the eigenvalue A, and w,. By examining this equation, one can see that the 
number of control variables (m) determines the dimension of the subspace in which the 
achievable eigenvectors must reside. 

Values of w, that yield an achievable eigenspace that is as close as possible in a least 
squares sense to a desired eigenspace can be determined by defining a cost function 
associated with the i^ 1 mode of the system 

^ = |(v.,-v 4 )"a,(v,-v Ji ) (2.23) 

for i = 1,...,/ where v a - is the i^ 1 achievable eigenvector associated with eigenvalue A, , is 
the i^ 1 desired eigenvector, and Qj i is an n-by-n symmetric positive semi-definite weighting 
matrix on eigenvector elements t . This cost function represents the error between the 
achievable eigenvector and the desired eigenvector weighted by the matrix Qd v 


t Superscript H denotes complex conjugate transpose (Strang 1980). 
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Values of w. that minimize 7, are determined by substituting (2.22) into the cost 
function for v a( , taking the gradient of 7, with respect to w, , setting this result equal to zero, 
and solving for w, . This yields 


w i = tfQd, L i I' 1 2 3 4 v rfi (2.24) 

where 

= [*■<,,!*- AY' & (2.25) 

and is the i th desired eigenvalue of the closed-loop system. Note in this development 
hij cannot belong to the spectrum of A. 

By concatenating the individual w,'s column-wise to form W and v aj .'s column-wise to 
form V a , equation (2.21) can be expressed in matrix form by 

W = [I m -GNY'GMV a (2.26) 

The feedback gain matrix that yields the desired closed-loop eigenvalues and achievable 
eigenvectors is given by 


G = W[MV U +AW]-' 


(2.27) 


Design Algorithm 

A feedback gain matrix that yields a desired closed-loop eigenspace is determined in the 
following way: 

1) Select desired eigenvalues , desired eigenvectors v ( j-, and desired eigenvector 
weighting matrices Qj i . 

2) Calculate w, 's using equation (2.24) and concatenate these column-wise to form W. 

3) Calculate achievable eigenvectors v a ■ 's using equation (2.22) and concatenate these 
column-wise to form V a . 

4) The feedback gain matrix G is then calculated using equation (2.27). 


Example 

An example will be presented to demonstrate the method. The design model is the 
lateral/directional dynamics of a high performance aircraft at low angle-of-attack. The 
model is based on a steady-state one g trim flight condition of forward cruise speed 
equaling 598 feet/second at 25,000 feet. It includes the four standard lateral-directional 
rigid-body degrees of freedom. The design goal is to improve the flying qualities by 



placing eigenvalues at level one locations and choosing eigenvectors to decouple the roll and 
dutch roll modes. 


The model is as follows: 

x = Ax + Bu (system dynamics) 
z = Mx + Nu (system measurements) 
u = u p + u c (total control input) 

u c — Gz (feedback control law) 
u p = G p S p (pilot control input) 


The system states are: 


x T =[fip,r s 


where 

P = sideslip angle (rad) 
p s = stability axis roll rate (rad/sec) 
r s = stability axis yaw rate (rad/sec) 
<f> = bank angle (rad) 


The system controls are: 


uT =l a roU a ya W f 


where 

a ro a — stability axis roll acceleration (rad/sec^) 
dyaw = stability axis yaw acceleration (rad/sec^) 

The measurements considered for feedback are: 

z T = [ P, r s d y $ ] T 


where 

p s = stability axis roll rate (rad/sec) 
r s = stability axis yaw rate (rad/sec) 
a y = lateral acceleration at the c.g.(g’s) 

P = sideslip rate (rad/sec) 


(2.28a) 

(2.28b) 

(2.28c) 

(2.28d) 

(2.28e) 

(2.29) 


(2.30) 


(2.31) 



The open-loop system matrices and open-loop eigenvalues and eigenvectors are given in 
Figure 1 . In this Figure, the open-loop eigenvectors have been scaled to allow comparison 
with the desired closed-loop eigenvectors. 

For this example, there are four states, two controls, and four measurements; therefore 
using DEA one can place four eigenvalues to desired locations and exactly place two 
elements of each associated eigenvector. The desired eigenvalues are chosen to yield good 
flying qualities. The desired roll and dutch roll eigenvectors are chosen to decouple the roll 
and dutch roll modes in the roll rate and sideslip responses. The desired eigenvalues and 
eigenvectors are given in Table 1 . In this Table, an x denotes eigenvector elements that are 
not weighted in the cost function. Therefore, the desired value for these elements is taken as 
arbitrary. Diagonal weighting matrices were used. Desired elements are weighted unity and 
other elements were weighted zero. 

The gain matrix to obtain the achievable eigenspace in the closed-loop system and 
closed-loop eigenspace is given in Figure 2. In this Figure, the closed-loop eigenvectors 
have been scaled to allow comparison with the desired eigenvectors. 

The closed loop eigenvalues have been placed at desired locations. As can be seen by 
examining the desired and achievable eigenvector elements, all the desired elements were 
obtained. 


DEA Conclusions 

DEA is a control synthesis technique for directly determining measurement feedback 
gains that will yield an achievable closed-loop eigenspace. For an observable controllable 
system that has n states, m controls, and / measurements one can determine a gain matrix 
that will place / eigenvalues to the desired locations and their associated eigenvectors as 
close as possible in a least squares sense to desired eigenvectors. 

Using DEA the designer has no direct control over augmentation gain magnitudes. 
Often it is not clear how to adjust the desired eigenspace in order to reduce individual 
undesirable gain magnitudes. To reduce undesirable gain magnitudes the designer must 
rely upon a strong physical insight into the dynamics or is forced to iterate on the design. 

The next section presents an eigenspace assignment methodology that overcomes this 
limitation. 


3.0 GAIN WEIGHTED EIGENSPACE ASSIGNMENT 

This section presents the development of the Gain Weighted Eigenspace Assignment 
(GWEA) methodology. This method allows a designer to place / eigenvalues at desired 
locations and trade-off the achievement of desired eigenvectors versus feedback gain 
magnitudes. 

The GWEA formulation builds upon a matrix formulation of DEA. The matrix DEA 
formulation is presented first followed by the GWEA formulation. The following 
development assumes complex matrices have been converted to real Jordan form (Broean 


Matrix DEA Formulation 

Given the observable, controllable system 

x = Ax + Bu (3.1a) 

where x € R" and u e R m , with system measurements available for feedback given by 


z = Mx + Nu (3.1b) 

where z e R( 

The total control input is the sum of the augmentation input u c and pilot's input u p 

u = u p + u c (3.2) 

The measurement feedback control law is 


u c = Gz (3.3) 

Solving for u as a function of the system states and pilot's input yields 

u = [I m - GNY' GMx+ [I m - GN]~ l u p (3.4) 

The system augmented with the control law is given by 

x = (A + B[I m - GN]~ l GM)x + B[I m - GNT' u p (3.5) 

The spectral decomposition of the closed-loop system is given by 

(A + B[I m - GNY' GM )v, = Vi (3.6) 

for i = 1 where A, is the i th system eigenvalue and v, is the associated i th system 
eigenvector. In matrix form this is given by 

(A + B(I m -GNy'GM)V = V A (3.7) 

Let W be defined by 


W = (I m -GNy'GMV 


(3.8) 


Substituting this equation into (3.7) and taking the vector value of each term yields 
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vec(AV)+vec(BW) = vec(VA) (3.9) 

where vec(X) is a vector valued function denoting a vector description of the matrix X 
(Graham 1981). Applying properties and rules for Kronecker products, equation (3.9) can 
be rewritten as 


(I,® A)vecV+ (I,® B)vecW = (A 7 ® I n )vecV (3.10) 

where ® denotes Kronecker product. Some definitions and rules for Kronecker products 
used in this development are given in the appendix. 

Solving for vecV as a function of A and vecW one obtains 

vec V = [(A r ® /„)-(/, ® A)] -1 (7, ® B)vecW (3. 1 1) 

or 


vecV=A D vecW 


(3.11b) 


where 


A d = [(A 7 ® /„)-(/, ® A)Y l (I, ® B) (3. 12) 

This equation (3.1 1) describes the achievable eigenvectors of the closed-loop system as 
a function of the desired closed-loop eigenvalues and vecW. If one could calculate a vecW 
that would make the achievable eigenvectors as close as possible in a least squares sense to 
some desired eigenvectors, it could be used to determine a gain matrix that would yield this 
eigenspace. One way this can be done is to define a cost function 


J e = vec T (V a - V d )Q d vec(V a - V d ) = e 7 Q d e (3. 13) 


where 


V a = achievable system eigenvectors 
V d = desired system eigenvectors 

Q d = symmetric positive semi-definite weighting matrix (eigenvector weighting matrix) 
e = vec(V a -V d ). 

This cost function represents the error between the achievable eigenvectors and the desired 
eigenvectors weighted by the matrix Qd- 

The value of vecW that will minimize J e can be obtained by taking the partial of J e with 
respect to vecW. This partial can be determined in the following way. Substituting equation 
(3.1 1) into equation (3. 13) for V a yields J e as a function of vecW. By applying the vector 
chain rule property (Graham 1981), the partial of J e with respect to vec W is given by 



dvecW 

Differentiating each term yields 




-r(e T Q d e) = Q d e + Q d e = IQ^e 
de 


(3.15) 


and 


de 


9 ( vec(V a - V „)) = —^—(A D vecW-vecV d ) = A T D (3.16) 


dvecW dvecW' a a " dvecW 
Therefore, the partial of J e with respect to vecW is given by 


dJ, 


dvecW 


hr. = 2 A T D Q d e = 2 A T D Q d (vec V a - vec V d ) 


(3.17) 


In this case, a closed form solution for vecW exists. This is obtained by setting equation 
(3.17) equal to zero and solving for vecW. This yields 


vecW = [A‘ D Q d A D Y'A T D Q d vecV d 


(3.18) 


where 


A d = [(A r d (8) /„)-(/, <8> A)] -1 (/, ® B ) 

and Ad is a block diagonal matrix of desired closed-loop eigenvalues. Note in this 
development the desired closed-loop eigenvalues cannot belong to the spectrum of A. 

The gain matrix that will yield the desired eigenvalues and achievable eigenvectors is 
obtained by solving equation (3.8) for G. 


G = W(MV a + NWy 1 


(3.19) 


By noting that 


V a =(vec r / ( ®/J(/ ( ® A 0 )(/ ;J ® WKI'tovecI,) (3.20) 

the feedback gains can be expressed as a function of the desired closed-loop eigenvalues 
and Wby 


G = W(MK 2 (7 ( i ® W)K 3 + NWy' (3.21) 


where 
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K 2 =(vec T I, »/„)(/, ®A D ) 
K 3 = (/, ® vec I, ) 


Gain Weighted Eigenspace Assignment Formulation 

The Gain Weighted Eigenspace Assignment formulation extends the Direct Eigenspace 
Assignment formulation to allow trading off eigenvector placement versus gain magnitudes , 
while still maintaining desired closed-loop eigenvalue locations. This is accomplished by 
appending a scalar measure of gain magnitude that is a function of vecW to the cost 
function given in equation (3. 13), determining partials with respect to vecW, and solving for 
the optimal solution by numerical iteration (Fletcher 1963). 

For this development, the scalar measure of gain magnitude is chosen to be a weighted 
sum of the squares of all the individual elements of the feedback gain matrix. To maYim i?* 
design flexibility, the gain magnitude term is formulated to allow weighting individual 
elements of the feedback gain matrix. A gain magnitude cost function that allows this can 
be formed in terms of the vector value of G. 

The vector value of Gas a function of vecW is obtained by taking the vector value of 
equation (3. 19) and noting that 


K = (vec T I, ® /„ )(/, ® vec V a ) 

(3.22) 

W = (vec T I, ® l m )(/, ® vec W) 

(3.23) 


thus yielding 


vecG = [(A/V a + AW) r ®/ m ]vecW 

= [(M(vec T I, ® /„)(/, ® A D vec W ) + N(vec T I, ® ® vec W)y T ® I m ]vec W 

(3.24) 


The gain magnitude cost function is given by 


J g = vec T (G)Q g vec(G) = g T Q g g (3.25) 


where 

Q g = symmetric positive semi-definite weighting matrix (gain weighting matrix) 
g = vec( G). 

This cost function represents the sum of the square of the individual feedback gains, each 
weighted by an element of the diagonal of the matrix Q g . The value of vecW that yields 
minimum gain magnitudes while achieving the desired closed-loop eigenvalues is 
determined by minimizing this cost function. 



Trade-offs between achievement of desired system eigenvectors and minimizing gain 
magnitudes can be made by forming the composite cost function 


J = p t J t +p g J g (3.26) 

where p e and p g are scalar positive cost function weights on the eigenvector placement error 
and gain magnitudes, respectively. 

Because eigenvectors can be scaled by an arbitrary constant, a unique solution to this 
cost function (3.26) does not exist when p e is zero ( or in practice when p e is small 
compared to p g ). To ensure a unique solution for all values of p e and p g it is necessary to 
constrain the eigenvectors to be unique. This can be accomplished by forcing one element 
of each eigenvector to be a specific reference value. To be consistent with the eigenvector 
error term J e in equation (3.26), this specific value is chosen to be an element of each 
desired eigenvector. This equality constraint can be expressed in the form of a penalty 
function (Bryson 1975) as 


J r = vec T (V a - V d )Q,vec(V a - V d ) = e T Q r e (3.27) 


where 

Q r = symmetric positive semi-definite weighting matrix (reference weighting matrix) 
e = vec( V a - V d ). 

This penalty function will be referred to as an eigenvector reference constraint. It represents 
the error between an element of each achievable eigenvector and the corresponding reference 
element of the desired eigenvector. The weighting matrix Q r is chosen to weight one 
element of each desired eigenvector. 

Appending this penalty function to the cost function (3.26) yields 


J = P e e T Q d e + p g g T Q g g + p r e T Q r e (3.28) 

where p e , p g , and p r are scalar positive cost function weights on the eigenvector placement 
error, gain magnitudes, and the eigenvector reference constraint; respectively. With this cost 
function, trade-offs between achievement of desired system eigenvectors and minimizing 
gain magnitudes can be made by choice of values of p g and p e . To ensure a unique 
solution the eigenvector reference constraint weighting p r should be chosen to be very large 
compared to the values of p g and p e . 


Partials of 7 with respect to vecW 


The value of vecW that will minimize J can be obtained by determining the gradient of J 
with respect to vecW. The partial of J with respect to vecW, is given by 


dJ 


= P< 


dJ. 


+P, 


dJ. 


dvecW ' dvecW * dvecW T dvecW 


+ Pr 


dJ. 


(3.29) 



The partial of J e with respect to vecW is given by equation (3.17). 

The partial of J g with respect to vecW can be determined in the following way. 
Applying the vector chain rule yields 


Jh_ 

dvecW 



Differentiating each term yields 


^(i r e,s)=e,s+2;*=2e,s 


and 


(3.30) 


(3.31) 


dg _ BvecG 
dvec W dvec W 


dG . dG . 

vec — — wee : 

dw u dw 2l 


■.vec 


dG 

dw. 


-iT 


ml , 


(3.32) 


where Wy is the (ij)^ 1 element of W . 

By applying the matrix chain rule, the partial of G with respect to wy is given by 


dG d 

AT - + W"T') 


f dw' 




(MK 2 (I i2 <8> W)K 3 + AW)-' + W] 


-^-(MK 2 (I i2 ® W)K , + AW)' 1 


The first partial with respect to wy in this equation is given by 


f 






(3.33) 


(3.34) 


where Uij is a matrix of order m-by-l which has unity in the (i j)* position and all other 
elements are zero. 


The second partial in equation (3.33) is given by 
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-^—(MK 2 (I 2 0 W)K 3 + AW)-' 


\ 

y 


~(MK 2 (I i2 0 + AW)- 


—~(mk 2 (i i2 0 w)a : 3 + aw) 

cAv.. ' 


(m/s : 2 (/ j2 0 w)/r 3 + aw) 


-1 


where 


(3.35) 


—~(MK 2 (I i2 0 W)tf 3 + AW) = MK 2 

ij 




*> + "Ir 

<7W 

y 


XI dw 

— MK 2 (1 2 ® -r — )/ l 3 + Af 

<7VV ; , 


dWy 


= MK 2 (I l2 ®U ij )K , + NU v 
Therefore, the partial of G with respect to Wy is 


(3.36) 


r)C 

— = Uy(MK 2 (I lt 0 W 3 + AW )- 1 


- W{MK 2 {I {2 0 W)K 2 + NWY\MK 2 (I i2 0 t/ iy ) A - 3 + NUyXMK^If 0 IV)AT 3 + AW ) -1 
= U 0 (MV a + AW ) -1 - G(MK 2 (I, 0 wct7 # ) + NU ij ){MV a + AW )' 1 
= (7 m - GN)U ij (MV a + NWy 1 - GMK 2 (I, 0 vecU tj )(MV a + AW )' 1 


(3.37) 

Taking the vector value of both sides of this equation yields 

dG 

vec—— - ((MV e + NW)~ T 0 (7 m - GN))vecUy - ((MV a + NW)~ T ®GMK 2 )vecV, ® vecU 0 ) 

OW ij * 


(3.38) 


Substituting equation (3.38) into equation (3.32) yields 



j^ = [(WV' + NW)-'®(l m -GN))lv'cU n ivect/,,1 ••• \vecUjf 

~ [{(MV a + NW)~ t <g> GMK 2 )[vec(I' <g> vecU n ): ■■■ \vec{I l ^vecU ml )^ 

= [vecU u :vecU 2i : . ••• \vecU ml } T ({MV a + NW)~' ®(I m -GN) T ) 
-[vec(I t ®vecU u )': ••• :vec(I, ®vec U ml )\ T ({MV a + AW)" 1 ®(GMK 2 f) 

(3.39) 

By noting that 


Li =[vecU u ':vecU 2l : ••• ivecf/^] 


(3.40) 


and defining 


U = [vec(l, ® vecU u ) : v£>c(/ ; ® vec U 2l ) : ■ ■ ■ : vec(I, ® vec £/„,)] (3.41) 

the partial of g with respect to vecW is 
dg 


~^~W = (W. + AW )' 1 ® ( 7 » “ GN) t ) - U T ((MV a + AW)" 1 ® ( GMK 2 ) t ) 


(3.42) 


Therefore, the partial of J g with respect to vecW is 
dJ 

= 2 [(< MV '. + ® (/. - GNf)-U T ((MV, + AW)’ 1 ® (GMAr,) r )j2,s 

(3.43) 

The partial of J r is determined in a manner similar to the partial of J e . The partial of J r 
with respect to vecW is given by 

dJ, 


Yh^W ~ 2A ° Q,e = 2A ° Qr(VeC V “ ~ V€C V “ } ( 3 - 44 ) 


Design Algorithm 

^ A feedback gain matrix can be calculated using the GWEA algorithm in the following 

1) Select desired eigenvalues , desired eigenvectors , and desired eigenvector 
weighting matrix Qj. 



2) Select gain weighting matrix Q g and reference weighting matrix Q r . 

3) Select cost function weights (see (3.28)) on eigenvector error p e and gain magnitude p g . 

4) Set cost function reference constraint weighting p r = 100.0*max( p e , p g ). 

5) Calculate VP that minimizes cost function (3.28) using analytic gradients in (3. 17), (3.43), 
and (3.44) along with numerical optimization techniques in Fletcher 1963. 

6) Achievable eigenvectors V a . can be calculated using equation (3.20). 

7) The feedback gain matrix G is then calculated using equation (3.24) or (3.19). 

8) Adjust cost function weights p g and p e as necessary to achieve desired design trade-offs. 


Example 

An example will be presented to demonstrate the method. This example uses the same 
model presented in Section 2; the lateral/directional dynamics of a high performance aircraft 
at low angle-of-attack. The open-loop system matrices and open-loop eigenvalues and 
eigenvectors are given in Figure 1 . 

The desired eigenvalues and eigenvectors are chosen to be the same as in the previous 
example. The desired eigenvalues are chosen to yield good flying qualities. The desired 
roll and dutch roll eigenvectors are chosen to decouple the roll and dutch roll modes in the 
roll rate and sideslip responses. The desired eigenvalues, desired eigenvectors, and 
eigenvector weighting matrix are given in Table 1. The gain weighting matrix Q g was set to 
identity. The desired eigenvector elements weighted in die reference weighting matrix Qr 
are given in Table 2. The reference constraint weighting p r was set to 100.0*max(p e , p g ). 

The algorithm was implemented in MATLAB t and executed on a SUN SPARC 10 tt 
workstation. Designs were determined for ten values of p g /p e . Solutions for each design 
were obtained in less than five minutes and within 50 iterations. Values of p g /p e ,J e , and 
J g for these designs are given in Table 3. A plot of J e versus the square root of J g is given 
in Figure 3. Gain matrices, closed-loop system matrices, and closed-loop eigenspaces for 
four values of p g / p e (0.0, 0. 1 , 100.0, 1 .0e05) are given in Figures 4 - 7. In these Figures, 
the closed-loop eigenvectors have been scaled to allow comparison with the desired 
eigenvectors. 

Cost function weighting p g /p e = 0.0 yields the DEA solution - four eigenvalues at 
desired locations and their associated eigenvectors are as close as possible in a least squares 
sense to the desired eigenvectors. The solution for this weighting is the same as the DEA 
solution presented in Figure 2. The RMS gain magnitude for this design is 5.73. 

Cost function weighting p g /p e = 0.1 yields four eigenvalues at desired locations and 
their associated eigenvectors close to the desired eigenvectors with a RMS gain magnitude 
of 4.08. This is a RMS gain magnitude reduction of approximately 28% compared to the 
DEA solution. 


t MATLAB is a registered trademark of The Math Works, Inc. 
ft SUN SPARC 10 is a registered trademark of Sun Microsystems, Inc. 



Cost function weighting p g /p e = 1.0e05 yields four eigenvalues at desired locations 
and eigenvectors that minimize feedback gain magnitudes (as defined by equation (3.25)). 
The RMS gain magnitude for this design is 2. 13. This is a RMS gain magnitude reduction 
of approximately 63% compared to the DEA solution. 


4.0 CONCLUDING REMARKS 

This report has presented the development of the Gain Weighted Eigenspace 
Assignment methodology. This provides a designer with a systematic methodology for 
trading off eigenvector placement versus gain magnitudes, while still maintainin g desired 
closed-loop eigenvalue locations. This was accomplished by forming a cost function 
composed of a scalar measure of error between desired and achievable eigenvectors and a 
scalar measure of gain magnitude, determining analytical expressions for the gradients, and 
solving for the optimal solution by numerical iteration. For this development the scalar 
measure of gain magnitude was chosen to be a weighted sum of the squares of all the 
individual elements of the feedback gain matrix. To achieve a solution it was necessary to 
constrain the system eigenvectors to be unique. This was accomplished by appending a 
penalty function to the cost function. 

An example was presented to demonstrate the method. In this example it was shown 
that cost function weighting p g / p e = 0.0 yielded the Direct Eigenspace Assignment 
solution - closed-loop eigenvalues at desired locations and their associated eigenvectors are 
as close as possible in a least squares sense to desired eigenvectors. As the cost function 
weighting p g / p e tended towards infinity the solution yielded closed-loop eigenvalues at 
desired locations and eigenvectors that minimized feedback gain magnitudes. Solutions 
yielding achievable eigenvectors close to the desired eigenvectors could be obtained with 
significant reductions in gain magnitude compared to the Direct Eigenspace Assignment 
solution. 
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APPENDIX 


This appendix presents some definitions and rules used in the development of the Gain 
Weighted Eigenspace Assignment methodology. For more information and proofs of these 
properties the reader is referred to Graham 1981. 


Definition of the VEC Operator 

Given a matrix A of order m-by-n, a vector valued function of matrix A, denoted by vecA, 
is defined by 


vecA = 


a n 
<h i 


where fly denotes the (i j)^ element of the matrix A. The vector vecA is of order mn-by-l 


Definition of the Kronecker Product 


Given a matrix A of order m-by-n and a matrix B of order r-by-p, the Kronecker 
product of the two matrices, denoted by A ® B, is defined as 


A®B = 


a u B 
°21 & 


a, 2 B 

U 22 & 


<hn & 


[a mX B a ml B - a^Bj 

where aij denotes the (i j) th element of the matrix A. The matrix A ® B is of order mr-by- 
np. 


Some Rules for Kronecker Products 

A ® (B + C) = (A ® B) + (A <8> C) 

(A + B) ® C = (A ® C) + (B ® C) 

(A ® B) (C ® D) = AC <8> BD (provided matrix dimensions are compatible) 
(A <8> B) T = A T ® B T 

(A <8> B)* 1 = A -1 ® B _1 (subject to existence of matrix inverses) 
vec(A Y B) = (B T ® A) vecY 



Open- loop System Matrices 


A = 


-0.1305 

0.0003 

-0.9978 

0.0537 

-11.1989 

-1.5271 

0.6757 

0 

2.9944 

0.1152 

-0.1529 

0 

0 

1.0000 

0 

0 

-0.0033 

-0.0319 



1.1179 

-0.1941 



0.0096 

1.3527 



0 

0 



0 

1.0000 

0 

0 

0 

0 

1.0000 

0 

-1.2576 

0.0535 

-0.0462 

0 

-0.1198 

0.0003 

-0.9992 

0.0538 

0 

0 



0 

0 



-0.0614 

-0.0669 



-0.0033 

-0.0319 




Open- loop Eigenvalue Matrix 


-0.2072 + 1 . 6584i 
0 
0 
0 


0 

-0.2072 - 1 . 6584i 
0 
0 


0 

0 

-1.4004 

0 


Open- loop Eigenvector Matrix 


1.0000 


1 

-3.6375 + 

3 . 7717i 

-3 

0.2110 - 

1 . 5598i 

0 

2.5092 + 

1 . 8800i 

2 


0000 -0.0147 
6375 - 3 . 7717i 1.0000 
2110 + 1 . 5598i -0.0569 
5092 - 1 . 8800i -0.7141 


0 

0 

0 

0.0043 


0.0026 

0.0043 

0.0535 

1.0000 


Figure 1 - Open-Loop System Matrices and Eigenspace. 



Gain Matrix and Gain Magnitude 
G = 

-0.3096 -0.8200 -5.1642 -0.6697 

0.0230 0.1138 -0.6396 2.1203 

ga inmag = 

5.73 


Closed-loop System Matrix 
A+B*inv(Im-G*N) *G*M = 


-0.1916 

-0.0519 

4.2491 

0 


0.0044 

-2.5051 

0.0552 

1.0000 


-0.9324 

0.0209 

-2.8083 

0 


0.0501 

-0.0136 

0.1509 

0 


Closed- loop Eigenvalue Matrix 

-1.50 + 1 . 50i 0 

0 -1.50 - 1 . 50i 

0 0 

0 0 


0 

0 

-2.50 

0 


0 

0 

0 

-0.0050 


Closed-loop Eigenvector Matrix 


1.0000 

-0.0225 + O.OOOOi 
1.4036 - 1 . 6084i 
0.0075 + 0 . 0075i 


1.0000 

-0.0225 - O.OOOOi 
1.4036 + 1 . 6084i 
0.0075 - 0 . 0075i 


- 0.0000 - 0.0000 

1.0000 -0.0050 

-0.0168 0.0537 

-0.4000 1.0000 


Figure 2 - DEA Design Closed-Loop System Matrices and Eigenspace. 






Gain Matrix 


G = 


-0.3096 -0.8200 -5.1642 -0.6697 

0.0230 0.1138 -0.6396 2.1203 

ga inmag = 

5.73 


Closed-Loop System Matrix 
A+B*inv ( Im-G*N) *G*M 


-0.1916 

-0.0519 

4.2491 

0 


0.0044 

-2.5051 

0.0552 

1.0000 


-0.9324 

0.0209 

-2.8083 

0 


0.0501 

-0.0136 

0.1509 

0 


Closed-Loop Eigenvalue Matrix 


-1.5000 + 1 . 5000i 
0 
0 
0 


0 

-1.5000 - 1 . 5000i 
0 
0 


0 

0 

-2.5000 

0 


Closed-Loop Eigenvector Matrix 


1.0000 


1 

-0.0225 + 

O.OOOOi 

-0 

1.4036 - 

1 . 6084i 

1 

0.0075 + 

0 . 0075i 

0 


0000 - 0.0000 

0225 - O.OOOOi 1.0000 

4036 + 1 . 6084i -0.0168 

0075 - 0 . 0075i -0.4000 


-0.0050 


- 0.0000 

-0.0050 

0.0537 

1.0000 


Figure 4 - GWEA Design p g /p e = 0.0. 
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Gain Matrix 
G = 

-0.5988 0.1195 -3.6945 -0.2224 

0.3561 -0.9005 -0.8321 0.9629 

gainmag = 

4.08 


Closed-Loop System Matrix 
A+B*inv(Im-G*N) *G*M 


-0.1927 

-4.2145 

4.7782 

0 


-0.0053 

-2.6323 

0.4762 

1.0000 


-0.9384 

1.0924 

-2.6800 

0 


0.0520 

-0.0081 

0.0720 

0 


Closed-Loop Eigenvalue Matrix 


-1.5000 + 1 . 5000i 
0 
0 
0 


0 

-1.5000 - 1 . 5000i 
0 
0 


0 
0 

-2.5000 

0 -0.0050 


Closed-Loop Eigenvector Matrix 


1.0000 

-1.5835 + 0 . 5681i 
1.4418 - 1 . 5829i 
0.7172 + 0 . 3384i 


1.0000 

-1.5835 - 0 . 5681i 
1.4418 + 1 . 5829i 
0.7172 - 0 . 3384i 


-0.1044 

1.0000 

-0.2844 

-0.4000 


0.0148 

-0.0050 

0.0525 

1.0000 


Figure 5 - GWEA Design p g /p e = 0.1. 
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Gain Matrix 
G = 

-1.0725 0.4415 -0.8696 -0.0573 

0.3095 -0.5105 -1.0925 1.2173 

ga inmag = 

2.25 


Closed-Loop System Matrix 
A+B*inv(Im-G*N) *G*M 

-0.1778 -0.0017 -0.9450 0.0515 

-10.0658 -2.8715 1.5224 -0.0120 

4.8368 0.3527 -2.4556 0.0918 

0 1.0000 0 0 


Closed-Loop Eigenvalue Matrix 

-1.5000 + 1 . 5000i 0 0 

0 -1.5000 - 1 . 5000i 0 

0 0 -2.5000 

0 0 0 -0.0050 


Closed-Loop Eigenvector Matrix 

1-0000 1.0000 -0.0637 0.0083 

-3.4556 + 2 . 0378i -3.4556 - 2.0378i 1.0000 -0.0050 

1.5052 - 1 . 5652i 1.5052 + 1.5652i -0.1801 0.0530 

1.8311 + 0.47261 1.8311 - 0.4726i -0.4000 1.0000 


Figure 6 - GWEA Design p g /p e = 1.0. 
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Gain Matrix 


G = 


-1.0238 

0.1539 


0.3554 

-0.3809 


-0.2791 

-0.9694 


0.1474 

1.4704 


gainmag = 
2.13 


Closed-Loop System Matrix 
A+B*inv(Im-G*N) *G*M 


-0.1662 

-11.0128 

4.4611 

0 


0.0023 

-2.7160 

0.1667 

1.0000 


-0.9402 

1.2526 

-2.6228 

0 


0.0511 

-0.0054 

0.1097 

0 


Closed-Loop Eigenvalue Matrix 


-1.5000 + 1 . 5000i 

0 

0 

0 

0 

-1.5000 - 1 . 5000i 

0 

0 

0 

0 

-2.5000 

0 

0 

0 

0 

-0.0050 


Closed-Loop 

> Eigenvector Matrix 



1.0000 


1.0000 


-0.0301 

0.0068 

-3.7635 + 

3 . 0195i 

-3.7635 - 

3 . 0195i 

1.0000 

-0.0050 

1.5322 - 

1 . 5744i 

1.5322 + 

1 . 5744i 

-0.0940 

0.0532 

2.2610 + 

0 . 2480i 

2.2610 - 

0 . 2480i 

-0.4000 

1.0000 


Figure 7 - GWEA Design p g /p e = 1.0e05. 



Desired Eigenvalues 
Aj =- 0.005 
A 2 = - 2.5 

h,4 = - 1.5 ± 1.5 j (o)n =2.1 (rad/sec), £ = 0.7) 


(spiral) 

(roll) 

(dutch roll) 


Desired Eigenvectors 


Mode 1 


(spiral) 


0 


Mode 2 


(roll) 


0 


Mode 3 
(dutch roll) 

1.0 


Mode 4 
(dutch roll) 

1.0 


1 x 0.0075 + 0.0075 j 0.0075 - 0.0075 j 

(x denotes eigenvector elements that are not weighted in the cost function) 


Table 1 - Desired Eigenvalues and Eigenvectors. 



Mode 1 

Mode 2 

Mode 3 

Mode 4 

(spiral) 

(roll) 

(dutch roll) 

(dutch roll) 

0 

0 

1 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 


Table 2 - Reference Weighting of Desired Eigenvector Elements. 



Pe 

Pg 

Pg ! Pe 

Sqrt( J g ) 

Je 

1.00 

0.0 

0.0 

5.7274 

1 . 99e-17 

1.00 

0.010 

0.01000 

5.3075 

0.00767 

1.00 

0.0667 

0.0667 

4.5356 

0.2926 

1.00 

0.10 

0.100 

4.0780 

0.6168 

1.00 

0.125 

0.125 

3.7757 

0.8785 

1.00 

0.20 

0.200 

3.1843 

1.5105 

1.00 

0.50 

0.500 

2.4804 

2.6575 

1.00 

1.00 

1.00 

2.2663 

3.3091 

1.00 

2.00 

2.00 

2.1821 

3.9256 

1.00 

10.0 

10.0 

2.1352 

4.7903 

0.010 

1.00 

100.0 

2.1319 

5.0980 

0.0010 

1.00 

1000.0 

2.1318 

5.1330 

1 . 00e-05 

1.00 

1 . 00e05 

2.1318 

5.1372 


Table 3 - GWEA Design Points . 
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